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Abstract 

LISA is an array of three spacecraft in an approximately equilateral triangle configuration which 
will be used as a low-frequency gravitational wave detector. We present here new generalizations of 
the Michelson- and Sagnac-type time-delay interferometry data combinations. These combinations 
cancel laser phase noise in the presence of different up and down propagation delays in each arm 
of the array, and slowly varying systematic motion of the spacecraft. The gravitational wave sen- 
sitivities of these generalized combinations are the same as previously computed for the stationary 
cases, although the combinations are now more complicated. We introduce a diagrammatic repre- 
sentation to illustrate that these combinations are actually synthesized equal- arm interferometers. 
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The Laser Interferometer Space Antenna (LISA) [1] is a space-borne gravitational wave 
(GW) detector mission which will use coherent laser beams exchanged between three widely 
separated spacecraft to study low-frequency (1CT 4 — 1 Hz) GWs. Modeling each spacecraft 
as carrying lasers, beam splitters, photo-detectors and drag-free proof masses on each of two 
optical benches, it has been shown that the six measured time series of Doppler shifts 

of the one-way laser beams between spacecraft pairs, and the six measured shifts between 
adjacent optical benches on each spacecraft, can be combined, with suitable time delays, to 
cancel the otherwise overwhelming phase noise of the lasers (Ais/v ~ 10~ 13 ) to a strain level 
h ~ 1CT 23 . This technique is called time-delay interferometry (TDI). 

Initial analyses of TDI □ □□□□ were for a non-rotating, rigid LISA array. However,the 
actual LISA orbits j?l [| produce annual rotation of the array. Recently, TDI applied to 
a rotating LISA was considered [9]. It was shown that the Sagnac effect leads to non- 
negligible light-time differences for light traveling around the array in the clockwise and 
counterclockwise senses. New candidate implementations of TDI involving both one-way 
measurements and laser phase-locking were presented which obviate the difficulties of a 
rotating array. In parallel work, Cornish and Hellings Q| pointed out that in addition to 
the rigid body rotation, more general inter-spacecraft velocities can occur (so called "flexing" 
of the LISA constellation). It was shown that this flexing will introduce phase noise at an 
unacceptable level, given the present laser frequency noise specification. 

In this paper we first generalize the original Michelson TDI combinations to an array 
with systematic spacecraft velocities, showing that these generalizations effectively cancel 
all laser phase noises (residuals second order in v/c, much smaller than the secondary LISA 
noise levels). We then analyze the generalized Sagnac combinations, showing that they too 
cancel laser phase noise adequately in a rotating and/or shearing LISA array. Finally, we 
emphasize that the gravitational wave sensitivities of these generalized combinations are the 
same as those for the non-rotating case. 

There are six beams exhanged between the LISA spacecraft, together with the six phase 
measurements Sij = 1,2,3) recorded when each transmitted beam is mixed with the 
laser light of the receiving optical bench. The phase fluctuations from the six lasers, which 
need to be canceled, can be represented by six random processes pij, where is the phase 
of the laser in spacecraft j on the optical bench facing spacecraft i. In what follows we 
assume the center frequencies of the lasers are all the same. 
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Since the LISA triangular array has systematic motions, the two one-way light times 
between any spacecraft pair are not the same Q . Delay times for travel between the space- 
craft must now be accounted for depending on the sense of light propagation along each 
link when combining these data. As before, we label the arms with single numbers given by 
the opposite spacecraft; e.g., arm 2 (or 2') is opposite spacecraft 2. We use primed delays 
to distinguish light-times taken in the counter-clockwise sense and unprimed delays for the 
clockwise light times. (Note that we have changed the data labeling conventions from that 
used in all previous papers by the last three authors. The subscript labeling of s^- in this 
paper is that of |9].) Explicitly: S23 is the one-way phase shift measured at spacecraft 3, 
coming from spacecraft 2, along arm 1. The laser phase noise in s 2 3 is P32{t — Li) — P23{t), 
where we take c = 1, so that L\ is the light time in the direction from spacecraft 2 to space- 
craft 3. Similarly, s^i is the phase shift measured on arrival at spacecraft 2 along arm 1' of 
a signal transmitted from spacecraft 3. The laser phase noise in s 32 is P23(^ — -^1) — ^32(^)5 
where L t is the light time in the sense from 3 to 2 along arm l'. Due to the relative mo- 
tion, L\ 7^ L 1 in general. For the further delays used in the TDI combinations we use the 
same conventions, being careful to distinguish light travel along arms with primes or not, 
depending on the sense of the measurement. For example, our notation for delaying the 
time series S32(t) by the clockwise light time in arm 1 would be 532,1 while delaying by the 
counterclockwise light time in arm l' would be 532,1'. As before, we denote six further data 
streams, Ty = 1,2,3), as the intra-spacecraft metrology data used to monitor the mo- 
tion of the two optical benches and the relative phase fluctuations of the two lasers on each 
of the three spacecraft. 

Cornish and Hellings have formulated TDI when the two delay times on each link, 
e.g. Li and L 1 are not only different (pure rotation) but also themselves functions of time. 
In the subscript notation for delays the order of the subscripts now becomes important for 
laser phase terms. The subscripts can no longer be permuted freely to show cancellation of 
laser noises in the TDI combinations and we will use a semicolon, instead of a comma, to 
emphasize this. (The other, secondary, noises in LISA are so much smaller, and the rotation 
and systematic velocities in LISA are so intrinsically small, that index permutation may still 
be done for them.) We will then go to first order expansions of the velocity, L, dropping 
quadratic terms in L/c and acceleration terms. This iterated time delay method, to first 
order in the velocity, is illustrated abstractly as follows. Given a function of time \l/ = ^(t), 
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time delay by Lj is denoted with the standard comma notation: 

V ti = *(t-Li(t)) (1) 

We then impose a second time delay Lj(t): 

= VQk-LjW-Litt-L^t)) 
~ ty(t — Lj(t) — Li(t) + Li(t)Lj) 

- + iftfLiLj (2) 

A third time delay Lk(t) gives: 

m, ijk = V(t - L k {t) - Lj(t - L k {t)) - Li(t - L k {t) - Lj(t - L k (t)))) 

~ y fc + yfcfL^Lj- + L fc ) + LjLfc] (3) 

and so on, recursively; each delay generates a correction proportional to its rate of change 
times the sum of all delays coming after it in the subscripts. 

Consider a constant-length equal-arm interferometer with one-way readouts for relative 
phase not only at the (central) spacecraft 1, but also at outlying spacecraft 2, and 3, with 
L 2 = L 2 = L 3 = L' 3 = L. The phase data s\ 2 , s 2 i, Si 3 , s 3 i can be combined to give the 
Michelson response 

S(t) = [s 31 + s 13)i ] - [s 2 i + s 12iL ] , (4) 

where we have purposely grouped the terms in square brackets to indicate that they provide 
the synthesized two-way phase data from each arm measured at the central spacecraft 
With these conditions, the laser phase noise is eliminated in S. An impulsive GW (duration 
short compared with L) will appear four times in the laser-noise-free time series S. 

With the unequal arm-lengths of space based interferometers, the time series S no longer 
cancels the laser phase fluctuations, leaving behind a remaining laser noise term that is 
proportional to the difference in the arm lengths to first order. In order to correct for 
this problem, it was shown in [2j that the two-way measurements entering into S must 
be differenced, with suitable delays, to again eliminate laser noise. This gives the TDI 
combination X Q: 

X = [(s 3 i + Si 3j2 ) + (S 2 1 + Si2,3'),22'] - [(S21 + S 12) 3>) + («31 + Si 3)2 ),33'] (5) 
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FIG. 1: Schematic diagram for X, showing that it is a synthesized zero-area Sagnac Interferometer. 
The optical path begins at an "x" and the measurement is made at an "o" 

together with the analogous time series, Y and Z, centered on spacecraft 2 and 3 obtained 
from equation (JHJ) by permutation of the indices. The response function of X to a gravi- 
tational wave pulse has twice as many terms as does S, and results in an "8 pulse" GW 
response. 

Equation (jSJ) shows that X is the difference of two sums of phase measurements, each 
corresponding to a specific light path (the continuous and dashed lines in Figure 1). The 
continuous line, corresponding to the first square-bracket term in equation (JHJ), represents a 
light-beam transmitted from spacecraft 1 and made to bounce once at spacecraft 3 and 2 
respectively. Since the other beam (dashed line) experiences the same overall delay as the 
first beam (although by bouncing off spacecraft 2 first and then spacecraft 3) when they are 
recombined they will cancel the laser phase fluctuations exactly, having both experienced 
the same total delays (assuming stationary spacecraft). For this reason the combination 
X can be regarded as a synthesized (via TDI) zero-area Sagnac interferometer, with each 
beam experiencing a delay equal to (L2 + L' 2 + L3 + L' 3 ). In reality there are only two beams 
in each arm (one in each direction) and the lines in Figure 1 represent the paths of phase 
information rather than paths of distinct light beams. 

If LISA has pure rotation, so that L 2 , L 2 , etc. are different but still time- independent, the 
original TDI combinations can be easily modified by being careful with the prime-noprime 
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notation. This has been done in the above equations for X. In this case these modified TDI 
combinations still cancel all laser noises exactly [lo| . 

In the general case of time- dependent spacecraft separations, in which systematic velocities 
enter as a perturbation on the static LISA unequal arm configuration, it has recently been 
pointed out by Cornish and Hellings n| that, for the LISA laser stability specifications of 
~ 30 Hz/ vTIz, the orbital variations, Li and bring in laser phase noise in the combination 
X at levels higher than the secondary noises in the lower part of the LISA frequency band 
(10~ 4 < / < 10~ 3 Hz). The same effect also prevents perfect cancellation of the phase noise 
in the other stationary-spacecraft TDI combinations . Improving the laser noise performance 
by about one order of magnitude, and reducing the LISA arm lengths were suggested as 
possible solutions to recover the required LISA sensitivity. Technical challenges may make 
further frequency stabilization down to this level unfeasible. In the following we show that 
there is another solution that requires no modification of the LISA hardware: it is possible 
to generalize TDI combinations (X, Y, Z) to remove the velocity dependence of the laser 
noise term remaining in the X combination. This is accomplished by further differencing 
the synthesized two-way measurements from each arm according to the diagram shown in 
Figure 2. As an example, let us assume that the velocity of spacecraft 3 relative to spacecraft 
1 and 2 is as shown in Figure 2. In this configuration the effective optical paths of the two 
synthesized beams (continuous and dashed lines) can be described as follows. One of the 
beams (continuous line, for instance) is first made to bounce off spacecraft 3 once, then 
spacecraft 2 twice, and finally makes one more bounce off spacecraft 3 before the phase 
measurement is made. Symmetrically, the other beam (dashed line) is first made to bounce 
off spacecraft 2 once, then spacecraft 3 twice, and once again off spacecraft 2 before the 
phase measurement is made. By delaying the beams in this manner we are able to average 
out the changes of the arm lengths taking place over a round-trip-light-time, making the 
two optical paths of the two beams essentially equal. 

The diagram shown in Figure 2 can be converted into a specific linear combination of the 
inter-spacecraft one-way phase measurements, s^, and metrology measurements performed 
onboard each spacecraft r^-. This new TDI combination, X±, is very insensitive to "flexing", 
and is given by: 
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FIG. 2: Schematic diagram for the new TDI combinations X\. Snapshots of spacecraft are shown 
at four times. Note that by properly delaying the beams within the two arms it is possible to 
average the effects of constant relative velocities of the spacecraft and equalize the optical paths 
of the two beams (continuous and dashed lines). 



X\ = [(s 3 i + Si 3;2 ) + (s 2 l + Si2;3');2'2 + (^21 + Sl2;3');33'2'2 + («31 + Sl3;2);33'33'2'2j 
~~ [( s 21 + Si2;3') + (^31 + Sl3;2);33' + ( s 31 + s 13;2)2'233' + ( s 21 + Sl2;3');2'22'233'J 
+ 2 [( T 21 ~~ T 3l) — ( r 21 — r 3l);33' ~~ ( r 21 — r 3l);2'2 + ( r 21 — r 3l) ;33'33'2'2 
+ (r 2 l — T3i) ;2 '22'233' — ( r 21 ~ ^31 ) ;2'233'33'2'2] (6) 

(with X 2 and X 3 obtained by cyclic permutation of the spacecraft indices.) Substituting 
into equation © the laser phase noise terms entering the and r^, and applying the 
expansion rules of equations (JH-BJ), it can be shown that, to first order in the systematic 
relative velocities of the spacecraft, laser phase noise is once again eliminated. This degree 
of suppression puts the laser noises several orders of magnitude below the secondary noises, 
even for a rotating/ shearing LISA array. Ultra-stable oscillator (USO) noises, which enter 
in the phase measurements of Xi, can be calibrated and removed in the same manner as for 
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FIG. 3: Schematic diagram for the generalized TDI combination ct\. By making each of the two 
beams circulate clockwise and counterclockwise once, it is possible to make their optical paths 
(continuous and dashed lines) essentially equal. (See text.) 

0. 



X 



13, in which 



Note that if one were to implement the locking configuration described in 
one of the lasers in spacecraft 1 is the master and the remaining five are slaved to it (so 
s i2 = S13 = r 21 — r 31 = r 32 — r 12 = r 13 — r 23 = 0), the expression above for X\ would reduce 
to 



Xil oc ked — [S31 + S21;2'2 + S21;33'2'2 + S 3 i ;33 ' 33 ' 2 ' 2 ] 

— [ s 21 + S3l ; 33' + S31;2'233' + S21;2'22'233'] > (7) 

where the two data sets s 3 i, s 2 i should now be regarded as two-way measurements. 

As with X±, the solution to "flexing" for the Sagnac combination a can be obtained by 
analyzing the diagram shown in Figure 3. Consider the sequence of one-way measurements 
starting from spacecraft 1, and let us propagate counterclockwise and clockwise the (syn- 
thesized) beams, represented by the continuous and dashed lines in Figure 3. Once again, 
for this demonstration we have assumed spacecraft 3 to have a non-zero velocity relative to 
the other two (stationary) spacecraft. The continuous and dashed lines correspond to the 
synthesized optical paths of the light transmitted by spacecraft 1. As Figure 3 suggests, 
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this way of combining the one-way measurements compensates for the arm length imbal- 
ance generated by the motion of spacecraft 3 during the round-trip time, making the delays 
experienced by the two synthesized beams essentially equal. By reference to Figure 3, the 
expression for the generalized Sagnac observable oti is: 

a l — [ s 31 + s 23;2 + Sl2;12 + s 21;312 + S32;3'312 + Sl3;l'3'312] 

— [S21 + S32;3' + S 13;1'3' + S31 ; 2'1'3' + S23;22'1'3' + Sl2;122'l'3'] 
+ 2 [( T 21 ~ r 3l) - (7"21 - T"3l);2'l'3'312 + (^32 ~ 7~12);3' 

+ (T32 — Tl2);12 ~ (^32 ~ Tl2);3'312 ~ (^32 ~ T12) ;122'1'3' 
+ ( r 13 — r 23);2 + (Tl3 ~ T 23 ) ; i'3' — (r 13 — T 2 3) ; 22'1'3' 

— (T23 — Tl3);l'3'312] (8) 

This expression coincides with that first derived in [^] except that now the order in which 
the delays are applied to the one-way measurements is important (as shown by the presence 
of a semicolon in equation (jHJ)). Expanding this expression using the results of equations (P 
we find that, unlike Xi, the residual phase noise of the laser(s) at spacecraft 1 remains 
to first order in the systematic velocities. The residual phase noise in a± is: 



P2i,i23i'2'3' [(^1 + L' 2 + L' 3 ) (Li + L 2 + L 3 ) - (Lx + L 2 + L 3 ) (L[ + L' 2 + L' s )} (9) 

Fortunately, although first order in the relative velocities, the residual is small, as it involves 
the difference of the clockwise and counterclockwise rates of change of the propagation delays 
on the same circuit. For LISA, the remaining laser phase noises in «j, % = 1,2, 3, are several 
orders of magnitude below the secondary noises. LISA's GW sensitivity in ct\ is essentially 
the same as for a. 

Since these generalized Michelson- and Sagnac-type TDI combinations involve twice as 
many terms as those entering into the original (stationary array) and modified (rotating 
array) TDI observables, the requirements in accuracies and precisions in physical quantities 
needed to synthesize the TDI combinations (such as arm length knowledge, synchronization 
of the onboard clocks, etc.) change. Our preliminary analysis indicates that the arm length 
and onboard clocks synchronization accuracies (needed to suppress laser noise to below sec- 
ondary noise sources in these new TDI combinations) will be more stringent than those 
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previously estimated for the original TDI combinations by a factor of about 2. Further anal- 
ysis is needed on these requirements and on practical issues concerning the implementation 
of TDI for LISA. 

This research was performed at the Jet Propulsion Laboratory, California Institute of 
Technology, under contract with the National Aeronautics and Space Administration. 
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